A semiclassical model of light mesons 



F. Brau and C. Semay 
Umversite de Mons-Hainaut, Place du Pare 20, B-7000 Mons, BELGIQUE 

B. Silvestre-Brac 
Institut des Seienees Nueleaires, Avenue des Martyrs 53, 
F-38026 Grenoble- Cedex, FRANCE 
(Dated: February 2, 2008) 

Abstract 

The dominantly orbital state description is applied to the study of light mesons. The effective 
Hamiltonian is characterized by a relativistic kinematics supplemented by the usual funnel potential 
with a mixed scalar and vector confinement. The influence of two different finite quark masses and 
potential parameters on Regge and vibrational trajectories is discussed. 
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I. INTRODUCTION 



Potential models, in particular semirelativistic ones, have been proved extremely success- 
ful for the description of mesons been developed in 
order to solve numerically the spinless Salpeter equation used in most works devoted to the 
study of light mesons. Variational calculations performed in well chosen bases , and the 
three-dimensional Fourier grid Hamiltonian method [4], for instance, can provide accurate 
numerical results. Nevertheless, it is interesting to obtain analytical results as well. The 
dominantly orbital state (DOS) description has been developed in this direction 0,0]. In 
this approach, the orbitally excited states are obtained as a classical result while the radially 
excited states can be treated semiclassically. It is then possible to obtain information about 
family of states characterized by a high orbital angular momentum. 

In Ref. 0], the DOS model has been applied to study light-light and heavy-light meson 
spectra with a confinement potential being a mixture of scalar and vector components. In 
the case of one or two quarks with vanishing masses, it was shown that linear Regge (orbital) 
and vibrational (radial) trajectories are obtained for an arbitrary scalar-vector mixture but 
that the ratio of radial to orbital energies is strongly dependent on the mixture. These results 



have been extended in Ref. 



by considering the presence of a quark and an antiquark with 
finite identical masses and the introduction of a Coulomb-like interaction plus a constant 
potential in addition to the confinement. In the limit of small masses and small strength 
of the short range potential, or in the limit of large angular momentum, it is shown that 
all calculations can be worked out analytically. One result is that slopes of orbital and 
radial trajectories depend only on the string tension and on the vector-scalar mixture in the 
confinement potential. Moreover, small finite quark masses do not alter significantly the 
linearity, specially in the case of dominant vector-type confining potential. As expected, the 
Coulomb-like potential has no effect on trajectories but its influence on the meson masses is 
determined in the approximation of the DOS description. Lastly, it is shown that the ratio 
of radial to orbital trajectory slopes depends on the scalar-vector confinement mixture only. 

In this work, we consider the case of a quark and an antiquark with finite but different 
masses. The potential taken into account is the same as the one considered in our previous 
work, that is to say a linear scalar-vector mixed confinement potential supplemented by 
a Coulomb-like interaction. An expression for the square mass of a meson is obtained in 
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the limit of small quark masses and small strength of the short range potential, or in the 
limit of large angular momentum. The coefficients of the square mass formula cannot be 
obtained analytically but they are computed numerically. Nevertheless a coverage of all 
possible situations is performed. The square mass formula is established and discussed in 
Sec. inland Sec. Illll respectivelv. Some concluding remarks are given in Sec. IIVI 



II. THE MODEL 



A. The dominantly orbital state description 

A detailed explanation of the technique of the DOS description is given in Ref. |6|. So we 
just recall the basic ideas and we only focus on differences between our work and previous 
ones. In all our formulas, we use the natural units h = c = 1. 

We consider a system composed of two particles with masses mi and m2 interacting via 
a scalar potential S{r) and a vector potential V{r) which depend only on the distance r 
between the particles. In the center-of-mass laboratory, the classical mass M of the system 
characterized by a total orbital angular momentum J is given by 



M(r,p„ J) = ^/p2 + + (^^ + aiSir)f + Jp2 + :^ + {m^ + a25(r))" + V{r), (1) 



where pr is the radial internal momentum. The parameters ai and a2 indicate how the 
scalar potential S{r) is shared among the two masses. These quantities must satisfy the 
following conditions to ensure a good nonrelativistic limit: 



«i + «2 = 1 and lim = 0. (2) 

mi— >oo 

A natural choice is to take 

m2 mi 
«! = and a2 = . (o) 

mi + m2 mi + m2 

We will use this prescription in this work. In the following, it is assumed that mi < m2. 
A parameter (3 = mi/m2 is introduced to measure the mass asymmetry, and the heaviest 
uark mass is denoted by m. Contrary to the philosophy adopted in our previous paper 
, it is necessary to work with the hamiltonian formalism and to give up the lagrangian 
formalism. 
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The idea of the DOS description is to make a classical approximation by considering 
uniquely the classical circular orbits, that is to say the lowest energy states with a given J. 
This state is defined by r = ro(J), and thus dr/dt = and Pr = 0. Let us denote Mq{J) = 
M{r = rQ,Pr = 0, J). In order to get the radial excitations, a harmonic approximation 
around a classical circular orbits is calculated. If the harmonic quantum energy is given by 
fi(J), then the square mass of the system with orbital excitation J and radial excitation 
?T,(0,1,. . . ) is given by (see Ref. 0) 

M^{J,n) = M^{J) + Mo{J)n{J){2n + 1). (4) 

The harmonic approximation is relevant only if it can be assumed that i7(J) ^ Mo(J). We 
will see in the next section that the last formula will be naturally obtained in the light meson 
case studied in this paper. 



B. Application to mesons 

As it is done in previous works, this formalism is applied to the study of mesons. The 
long range part of the interaction between a quark and an antiquark is dominated by the 
confinement which is assumed to be a linear function of r. As its Lorentz structure is not 
yet determined, we assume that the confinement is partly scalar and partly vector. The 
importance of each part is fixed by a mixing parameter / whose value is for a pure vector 
and 1 for a pure scalar. The short range part of the interaction is assumed to be of vector- 
type and given by the usual coulomb-like potential. Thus we have 

Sir) = far, (5) 

in which a is the usual string tension, whose value should be around 0.2 GeV^, and 

V{r) = {l-f)ar-- (6) 

r 

in which k is proportional to the strong coupling constant a^. A reasonable value of k should 
be in the range 0.1 to 0.6. The quark masses appearing in a spinless Salpeter equation are the 
constituent masses. Their values are model dependent, but generally we have m„ = 171^ ^ 0.2 
GeV and ^ 0.5 GeV. 

We are interested only in Regge trajectories which give the behavior of in term of J 
for large value of J. So to compute the square mass formula for the particular interaction 
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we will use the same technique as in our previous study. A detailed explanation is 
given in Ref . Q . For a given set of parameters / and (3, we have to compute the value of tq 
as a function of J, in the limit of great values of J. In the general case, this radius is given 
by a polynomial equation of 20th degree. This equation is by no mean obvious to obtain 
and cannot be considered as a simple extension of the two equal mass case. The software 
Mathematica has been used to perform the longest calculations. The general formulas are 
very complicated and too long to be given in detail here, but they can be obtained on 
request. 

When /? = 1, the old formalism must be recovered. This is not trivial to demonstrate, but 
we have verified this point. In particular, when the two particles have the same mass, the 
degree of the general polynomial is reduced and the solution is one of the root of a second 
degree equation. In this case, the determination of ro and all subsequent calculations can 
be performed analytically and are identical of the ones performed in Ref. Q. 

When (5 is different from 1, a numerical solution for ro must be computed numerically 
for each value of parameters / and (3. Note that when / = 0, only one relevant solution (ro 
real positive) exists (see Ref. To find the physical relevant value of ro, we generate the 
solutions for increasing values of / from the known solution at / = 0, in order to be sure 
that the solution is a continuous function of /. This procedure is repeated for each value of 
(3. Once the radius is calculated, its value is used to compute the square meson mass. 

Finally, we obtain the square mass of the meson under the following form 



where 0</<l,0</9<l, and m is the mass of the heaviest quark. The values of 
coefficients A and as a function of the mixing parameter / and the asymmetry parameter 
(3 are given in Figs. ^ All coefficients are discussed in the next section. It is worth noting 
that the term in Vl^ naturally does not appear as it is of higher order. 

When /3 = 1, that is to say rrii = m2 = m, formula (jZj) reduces to formula (20) of Ref. Q], 
and coefficients A, B, C, D and E are given by analytical expressions which can be found 
in this reference. We checked theoretically as well as numerically that our new expressions 
tend towards the previous ones at the hmit /3 — 1. 




(7) 
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III. DISCUSSION OF THE MODEL 



Equation ((Tj) is the square mass of a meson for large value of J. Though coefficients 
of the formula cannot be obtained analytically, they are calculated for < / < 1 and for 
< f3 < 1, which covers all possible physical situations provided the mass m appearing in 
the formula is the mass of the heaviest quark. One can make the following comments. 

For large value of J, the dominant term of the square mass formula is linear in J. This 
explains the existence of linear Regge trajectories. The slope of the Regge trajectories 
depend on three parameters: the string tension a, the mixing parameter / and the mass 
asymmetry parameter f3. The string tension is a constant and gives the energy scale of 
the slope. The slope depends strongly on the value of /, which is assumed to be also 
an universal quantity independent of the system. Unfortunately, as we have seen in our 
previous work experimental data from symmetrical mesons (p and families) cannot 
provide reliable information about the value of /. The calculations favor slightly a dominant 
vector- type confining interaction (/ = 0). But this work shows that the slope depends also 
on the parameter j3 whose values is fixed by the system. On Fig. ^ one can see that the 
/3-dependence of the coefficient A(f, (3) is only marked when / > 0.5 and /3 < 0.5. It is 
also shown in our previous work [7| that the conclusions we can obtain from our square 
mass formula can be applied to light mesons only, that is to say mesons containing u, d ot s 
quarks. In this case, the only physical relevant value for f3 is given by the ratio mu/ms- For 
constituent masses, the value of this ratio must be around 0.5. This shows that the slope of 
the Regge trajectory for K* mesons cannot be sensitively different from the ones for p and 
(p mesons. This is specially verified if the confinement potential is dominantly of vector-type 
(/ ~ 0) as it seems favored from experimental data. 

One can see on Fig.Q]that the slope of the vibrational trajectories proportional to E{f, (3) 
presents the same characteristics as the coefficient A{f,(3). We can conclude that the slope 
of the vibrational trajectory for K* mesons cannot be sensitively different from the ones 
for p and mesons. Obviously, the ratio R{f,P) = 2 E{f, P)/A{f, P) is characterized by a 
similar behavior as a function of / and f3. 

A term proportional to ^/j deforms the Regge trajectories for low value of J. This term 
vanishes if m = 0, that is to say the two quarks are massless, or if /5 = 0, that is to say the 
lightest quark is massless. So, the linearity of the Regge trajectory is perfect if one quark is 
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massless. This is also the case if the confinement potential is a pure vector-type interaction 
(/ = 0) whatever the masses are. The value of B{f, /3) decreases rapidly with / and P value 
parameters. This fact implies that the linearity of the Regge trajectory is only broken in 
the case of dominantly scalar-type confinement interaction. As it is mentioned above, this 
situation is not the one favored by experimental data. Note that the term proportional to 
y/j is also independent of the strong coupling constant k. 

The position of the trajectories with respect to the energy axe stems from three contribu- 
tions; one C(/, /3)m^ is due to the finite mass the quarks, another negative one a D{f, (5)k is 
due to the strong coupling constant and the last one aE{f,P) refiects the zero point motion 
of the harmonic vibration. These terms arc of minor importance since the zero point energy 
of the orbital motion cannot obviously be calculated in our model. Nevertheless it is inter- 
esting to discuss a little bit the behavior of coefficients C and D. The coefficient C{f,f3) 
depends clearly on /3 whatever the value of /. When f — 0, this coefficient can be calculated 
exactly: B(0, P) — 4:{1 + 0^). This term is linked to the contribution of the quark masses 
to the meson total mass. The presence of the Coulomb-like potential decreases the mass of 
the meson but has no infiuence on the Regge trajectories, as it is expected for a potential 
which becomes very weak with respect to the confinement interaction as quark intcrdistance 
increases. Again the coefficient £>(/,/?) is only affected by the mass asymmetry for / > 0.5 
and (3 < 0.5. 

Finally, let us remark that there is no coupling between orbital and radial motion for 
large J values (absence of terms nJ). This is only a consequence of the Coulomb-|-hnear 
nature of the quark- antiquark potential. This may not be true for other types of potentials. 

IV. CONCLUDING REMARKS 

We have shown that all light mesons exhibit linear orbital and radial trajectories in the 
dominantly orbital states (DOS) model. Slopes of both type of trajectories depend only 
on the string tension, the vector-scalar mixture in the confinement potential and the mass 
asymmetry between the quark and the antiquark. 

In the limit of small quark masses, it turns out that the mass asymmetry between the 
quark and the antiquark alters sensitively the slope of the trajectories only for large mass 
asymmetry and for dominant scalar-type confining potential. As the experimental data 
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favor a vector-type confining potential (see our previous paper [7|) and as tlie mass ratio 
of non-strange quark over the strange quark is not very small, tlie slope of the Regge and 
vibrational trajectory for K* mesons cannot be sensitively different from the ones for p and 
mesons. This is in agreement with experimental meson spectra. 
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FIG. 1: Coefficient A{f, (3) and E{f, (3) of formula ((Tj) as a function of the mixing parameter / and 
the mass asymmetry parameter (5. 
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